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We present a simple scheme to determine the diffusion properties of a thin slab of strongly scattering 
material by measuring the speckle contrast resulting from the transmission of a femtosecond pulse with 
controlled bandwidth. In contrast with previous methods, our scheme does not require time measurements nor 
interferometry. It is well adapted to the characterization of samples for pulse shaping, non-linear excitation 
through scattering media and biological imaging. © 2013 Optical Society of America 



The propagation of optical waves in scattering media 
is both a fundamental and a very applied topicwhich re- 
quires a deep understanding of how a wave propagates 
in a scattering medium An important transport pa- 
rameter that governs it is the transport mean free path, 
It . Upon multiple scattering, that is when the dimension 
L of the system is much larger than tt, virtually no bal- 
listic light transport is present and for k£t ^ ^ (where k 
is the wavenumber) we can treat, in first approximation, 
the propagation as a diffusive process governed by the 
diffusion constant, D. 

Coherence plays an essential role in propagation. In 
the continuous-wave domain, characterized by a single 
pulsation loq, multiple scattering gives rise to a very uni- 
versal random interference figure: the speckle For a 
polarized monochromatic wave, the spatial distribution 
of the speckle intensity / follows the Rayleigh distribu- 
tion P{I) = exp(— //(/))/(/), where (/) is the mean 
intensity, and the contrast is unity. 

In a slab geometry, measuring scattering properties, 
in particular it, D or the energy velocity ve is by no 
means trivial. The angular width of the coherent back- 
scattering cone [sj , or, in the case of thin slab geometry, 
total transmission studies Q, give an evaluation of it- 
The spread of a short pulse speckle correlations in 
time [6] or frequency phase sensitive interferometric 
measurements [t] or variation of the effective refractive 
index ^ allow to extract the diffusion constant D. Si- 
multaneous measurement of D and it [9 10 give access 
to Ve. 

The relevant time, in the dynamical case, is the Thou- 
less time Td 11 , with associated bandwidth Avd ~ l/Td- 
The characteristic time of diffusion ^ L'^ /(tt'^D) is the 
time it takes for a photon to leave the medium, related 
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which 



to the diffuse traversal time r; 

takes also into account the direction of the propagation 
of light (r/ > Td). In essence, a short laser pulse of du- 
ration Tp < Ti transmitted through the medium will be 
scattered all over the medium and its duration will be 



extended by approximately ti by multiple scattering. 

Earlier works have studied how the speckle depends 
on the coherence of the source [Ts] or is modified when 
a monochromatic laser is frequency-tuned in the GHz 
range 14 . The temporal behavior of the speckle has been 



measured with a streak camera in the sub-ns range 15 



or via non-linear interferometry in the fs range 16 . Some 



related works on temporal aspects of speckle from surface 
scattering have also been studied (see for instance [l7]). 
They rely either on temporal measurement (direct or 
via non-linear interferometry or Fourier Transform) or 
by time integrated contrast measurement on a camera. 

In this paper, we introduce a simple method to deter- 
mine the diffusion properties of a thin slab of strongly 
scattering material, in a regime where the other tech- 
niques are not easily applicable. The method is based on 
a simple contrast measurement on a CCD camera and il- 
lumination with a femtosecond laser of variable spectral 
bandwidth. 

The investigated samples are 3D scattering media with 
slab geometry made by sedimentation through verti- 
cal deposition from a 5 wt% (weight- weight percentage) 
water suspension of ZnO powder as described in refer- 
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. This technique allows to grow 
nous samples with a very flat surface, and filling fraction 
around 50%. The particle size distribution is 230±70 nm. 

In the simple model of reference [2], the contrast of 
the speckle image is shown to be C = 1/a/]V, where N is 
the number of independent speckle patterns added inco- 
herently. In the present case, the medium is illuminated 
by a pulse with a bandwidth Avp larger than the band- 
width of the medium Ai>m ~ Avd- After propagation 
through the medium, there are N ~ Avp/Avm inde- 
pendent spectral speckles patterns, making the contrast 
C ^ 1/ y^Ti X Afp. This implies that a measurement of 
the contrast as a function of Ai^p will give access to r/, 
and therefore D. 

The experimental setup is shown in Figure [T] The 
light source is a mode-locked Ti:Sapphire laser (Spec- 
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Fig. 1. Experimental setup with variable slit to control 
the spectral bandwidth, typical measured spectrum and 
speckle image on the CCD camera 



tra Physics Mai Tai 90 fs pulse and 80 MHz repetition 
rate), with a central wavelength of 748 nm. The 2 mm 
coUimated TEMqo output of the laser is sent in a zero dis- 



persion line 19 in order to select the spectral bandwidth 



of the light before the scattering medium. It consists of 
a series of 3 optical elements, a grating (1200 pitch per 
mm), a variable slit attached to a mirror, and a lens be- 
tween the mirror and the grating at the focal distance 
/ = 50 cm from both elements. After the zero dispersion 
line, the beam is focused at the surface of the scatte- 
ring medium with a / = 50 cm lens. Transmitted light 
is collected with a microscope objective (Olympus lOOX, 
NA = 0.8) and the back focal plane of the objective is 
imaged with a / = 20 cm lens on a CCD camera (AVT 
PIKE F-IOOB). The object plane of the microscope is 
moved by ^ 10 /im from the sample surface in order to 
increase the size of the speckle grain on the CCD Cam- 
era. A linear polarizer is placed between the objective 
and the lens to ensure that a single polarization com- 
ponent is imaged. To monitor the input spectrum the 
laser is partially deflected onto a spectrometer (Ocean 
Optics HR4000, resolution < O.lnm) by a beamsplitter 
after the zero dispersion line. Spectrum and image are 
acquired simultaneously. The pulse bandwidth at the en- 
trance of the medium is between the bandwith of the 
laser AAmaa; = 8 nm and AA„ii„ ~ 0.3 nm, limited by 
diffraction on the slit. 

A typical image on the CCD is shown in inset on Fig.[T] 
The contrast of the speckle is calculated within a statis- 
tically homogeneous region of the CCD, either in the 
speckle image (speckle), either not (dark). The contrast 
C = CT (speckle) /((speckle) — (dark)) , where () and cr() 
are the mean and standard deviation of the pixel inten- 
sities within the defined region. Error of C is evaluated 
based on the standard deviation of the dark image. 

Samples with different thicknesses were illuminated 
with light pulses of various bandwidths lS.Vp in order 
to investigate the link between the speckle contrast and 
D. To quantify the contrast modification with Ai/p in a 



more complete formalism, we treat a thin strongly scat- 
tering sample using a slab geometry and assuming a non- 
absorbing medium. This treatment, however, is general 
and can be extended to any geometry and medium in- 
cluding those in which both scattering and absorption 
are significant. The first question is to determine the 
distribution of transit time through such a medium. A 
general treatment of this problem has been introduced 
in [20]. In the multiple scattering regime, the spatially 
integrated transmittance (t) as a function of time for 
negligible absorption reads (from Eq. (15) of 20 ): 



-{L + lt)e + (3i - lt)e ^ 



(3i + £j)^ ~1 



(1) 



We neglected contributions by more than 3 reflections as 
they do not contribute significantly to Tj^ in the numer- 
ical calculation. This time-dependent transmission has 
to be compared to the pulse duration in order to deter- 



mine the contrast. In reference 13 a formalism is intro- 



duced for this purpose, in the context of continuous- wave 
laser propagating through a very thick medium, which is 
also valid in the case of a short pulse. A spatially coher- 
ent source of spectrum S{X) going through a scattering 
medium of time distribution fL{t) = T[t)/ T{t)dt 
produces a speckle with a contrast C given by the fol- 
lowing integral: 
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and where / is, with a change of variable, a Fourier trans- 
form of TL{t), and reads: 
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dt. (3) 



It is important to note that all parameters of the problem 
are not fully independent as Ti(t) and C depend in a 
non-trivial way on all the parameters L, £j and D. As 
a consequence it is not possible to infer all parameters 
independently from contrast measurements. 

By changing the width of the variable slit, measure- 
ments of the contrast were made over a range of spec- 
tral bandwidths IS.Up. As expected (Fig. [2]). As expected 
for large values of the bandwidth the contrast asymp- 
totically decreases with increasing bandwidth in the 
form C 
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lS.Vp and decreases with Az/p faster for 



thicker material as the traveling time gets longer. The 
monochromatic limit where C = 1 is reached when 
Ai^p < l^Vdi i-G. when Tp is identical to the confinement 
time of the medium 21 . 

Once known C(Azyp), the diffusion constant D can be 
obtained if the transport mean free path It is measured. 
it can be extracted from a total transmission measure- 
ment by employing the optical analogue of the Ohm's 
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Fig. 2. Contrast of the speckle for 3 different thick- 
nesses of ZnO. Dots are the experimental measurements, 
black line is a calculation following the model developed 

by 



20 . Inset: Inverse of the transmission for ZnO layers 



of different thickness L in fj,m. 

law: the angular integrated total transmission T is in- 
versely proportional to the thickness L of the sample. A 
simple expression relating T to L is l22]: 



T{L) = 



Ig sinh(22:e/4)sinh(ze/4) 
Ze sinh((L + 2ze)/4) 



(4) 



where we have assumed symmetric boundary conditions 
and = ^ln(j±|^), zq = ^f^, 4 is the ab- 
sorption length and i? is a correction factor due to the 
diffuse reflectivity of the boundary conditions 23 , which 
we calculated to be i?=0.49 for a ZnO sample in air. The 
interface between the sample and the glass substrate is 
here neglected as index matching is almost achieved be- 
tween glass (n=1.46) and the sample ((n) ^1.5). Eq. [4] 
holds in the multiple scattering regime {it ^ L) which 
is the case of the samples under study. When absorption 
can be neglected, Eq. [i] simphfies to T{L) ^ ^ it j^- 

Total transmission as a function of sample thickness 
is measured by focusing a laser beam at 748 nm on a 
set of samples placed at the entrance of an integrating 
sphere. The integrated light transmitted at all angles is 
sent via an optical fibre to a spectrometer, which mea- 
sures the total transmission 10 . The inset in figure [2] 
shows a plot of the inverse of the transmission T{L) as a 
function of the thickness for A — 748 nm. Vertical error 
bars are based on the variation of transmission for differ- 
ent position on the sample while the horizontal ones are 
the standard deviation of the sample thickness upon re- 
peated measurements. Fitting this data set to the model 
gave it = 2.1 ± 0.2/im and £a > 100 /im such that ab- 
sorption is negligible. 

C^Avp), L and it being measured independently, by 
fitting Eq.[2]to all data, D is found to be {29 ± 'i)m'^ 
for our ZnO samples. Fits are shown in Fig. [2] for 3 sam- 
ples. This allows ve to be calculated for our medium as 
(0.47 ±0.04) X 10^ms"\ i.e. ve ^ 0.16 c, where c is the 
speed of light in vacuum. 

Alternatively, the description of the speckle contrast 
in terms of interference of N independent spectral modes 
provides a straightforward estimate of the diffusion prop- 



erties of the medium: using C ^ l/-y/r/ x Ai>p, we esti- 
mate the diffusion coefficient to be (19±2) m^s^^, close 
to the rigorous value previously calculated. 

In conclusions, using a CCD camera and a femtosec- 
ond light pulse, we demonstrate a simple set-up to char- 
acterize the diffusion properties of a scattering medium. 
In order to avoid complex interferometry measurement 
of time and phase we exploit the ability of a multiple 
scattering medium to mix spatial and spectral modes 
at the output 
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In view of the recently realized 
spatio-temporal control of light pulse in scattering me- 
dia 24p5 , our concept applies for a wide range of Thou- 
less time, and therefore allows to characterize the diffu- 
sion coefficient D for various scattering media, for pho- 
tonics and biological applications. 
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